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Abstract. In this talk I have presented the data analysis results of 
I I . extracting properties of halo WIMPs: the mass and the (ratios between 

' the) spin-independent and spin-dependent couplings/cross sections on 

nucleons by the AMIDAS website. Although non-standard astronomical 
^ ^1 setup has been used to generate pseudodata sets for our analyses, it has 

^ ' been found that, without prior information/assumption about the local 

, density and velocity distribution of halo Dark Matter, these WIMP 

' ' ■ properties have been reconstructed with ~ 5% to ^ 40% deviations 

, ' from the input values. 

> ■ 

CN 1 Introduction 

T— H 

In order to extract properties of halo WIMPs (Weakly Interacting Massive Par- 
Q^ ■ tides) by using data from direct Dark Matter detection experiments as model- 
I independently as possible, we have developed a series of data analysis method for 
reconstructing the one-dimensional WIMP velocity distribution function (Drees 
k Shan 2007) as weh as determining the WIMP mass (Drees & Shan 2008), the 
spin-independent (SI) WIMP coupling on nucleons (Shan 2011) and the ratios 
k>( ■ between different WIMP couplings/cross sections (Shan 2011). Moreover, in col- 
j_] ■ laboration with the DAMNED (DArk Matter Network Exclusion Diagram) Dark 
_C3 ; Matter online tool (DAMNED), part of the ILIAS Project (ILIAS), the "AMIDAS" 
(A Model-Independent Data Analysis System) website for online simulation/data 
analysis has also been established (AMIDAS; Shan 2010, 2009). 

In this article, in order to demonstrate the usefulness and powerfulness as 
well as the model-independence of the AMIDAS package for direct Dark Matter 
detection experiments, I will analyze blindly some pseudodata sets generated for 
different detector materials and present the reconstructed WIMP properties. This 
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means that I will simply upload these data sets onto the AMIDAS website and 
follow the instructions to reconstruct different WIMP properties without using any 
information about the input setup used for generating the pseudodata. For cases 
in which some information about WIMPs (e.g., the mass m^) and/or Galactic 
halo (e.g., the local Dark Matter density po) is required, I will naively use the 
commonly used/favorite values for the data analyses. 

After that I show the blindly rciconstructed properties of halo WIMPs in Sec. 2, 
in Sec. 3 I will reveal the input setup used for generating the analyzed data and 
compare the reconstructed results to them. Finally, I conclude in Sec. 4. 



2 Reconstructed WIMP properties 

In this section, I present the reconstructed WIMP properties analyzed by the 
AMIDAS website. While in each uploaded file there are exactly 50 data sets, in 
each data set there are on average 50 recorded events (i.e., 50 measured recoil 
energies)^; the exact number of total events is Poisson distributed. For simplicity, 
the experimental minimal and maximal cut-off energies have been set as and 
100 keV for all data sets. 

In order to check the effect of using a "wrong" elastic nuclear form factor, two 
forms have been considered for the SI WIMP-nucleus cross section in our analyses. 
One is the simple exponential form: 

F^iQ) = e-«/«° . (2.1) 

Here Q is the recoil energy transferred from the incident WIMP to the target 
nucleus, Qo is the nuclear coherence energy given by Qq = 1.5/mN-Ro' where 
Rq = [0.3 + 0.91 (mN/GeV)^^^] fm is the radius of the nucleus and mN is the 
mass of the target nucleus. Meanwhile, we used also a more realistic analytic form 
for the elastic nuclear form factor: 



3ji(gi?i) 



qRi 



(2.2) 



Here ji{x) is a spherical Bessel function, q = \/2rrmQ is the transferred 3- 
momentum, for the effective nuclear radius we use Ri = yJR\ - 5s2 with Ra ^ 
1.2 A^/^ fm and a nuclear skin thickness s ~ 1 fm. For the SD WIMP-nucleus 
cross section, we only used the "thin-shell" nuclear form factor: 

p. J Jo , for qR^ < 2.55 or qR^ > 4.5 , 

const. ~ 0.047, for 2.55 < qRx < 4.5 . 
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Results 

Reconstructed WEVIP mass 
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The reconsbucled WIMP mass and the up]>er and lower bounds of its \c slatisiical uncertainly, m^.n and m^.irliavB 
been eslimaled by Eqs. (34) and (40) of Ref. [2], respectively; while %,combineil by tbe j;^-litting defined in Eq. 
(51) of Ref. [2], whicb combines the estnnators for and m^^ with each other. The teconslructed WIM^g 

ombmed as well OS m^.n and m^.a shown here have been corrected by the algoridimic Qmax matchirrg & 
faiRe£[2]. 



Fig. 1. The WIMP mass reconstructed with a target combination of ^*Si + ^''Ge 
nuclei. Two forms of the elastic nuclear form factor given in Eqs. (2.1) and (2.2) have 
been used in the upper and lower frames, respectively. 



2.1 WIMP mass 

As one of the most important properties of halo WIMPs as well as the basic 
information for reconstructing other quantities in our model-independent analysis 
methods, I consider at first the determination of the WIMP mass by means of 
the method introduced in Drees & Shan (2008). 

In Figs. 1 I show the reconstructed WIMP masses and the upper and lower 
bounds of their la statistical uncertainties. The usual target combination of 
^^Si + ^^Ge nuclei has been used for this reconstruction, whereas two forms of 
the elastic nuclear form factor given in Eqs. (2.1) and (2.2) have been used for 
determining in the upper and lower frames, respectively. While m^^n with 
n = —1, 1, 2 and m^^cr have been estimated by Eqs. (34) and (40) of Drees & Shan 
(2008), respectively, m^, combined bas been estimated by the x^~fitting defined in 
Eq. (51) of Drees & Shan (2008), which combines the estimators for m^," ^^'^ "^x.o- 
with each other. The reconstructed WIMP mass combined as well as m^^n and 
m^.g. shown here have been corrected by the iterative Qmax^niatching procedure 



^Note that we considered here only data sets with pure WIMP signals, possible unrejected 
background events are neglected. 



4 



Title : will be set by the publisher 



described in Drccs & Shan (2008). 

It can be found here that, although all single estimators {m^^n with n = —1, 1, 
2 and m^^a) give generally a (relatively lighter) WIMP mass of ~ 50 GeV or even 
lighter and a la upper bound of ^ 130 GeV, the m,ean vahies of the combined (in 
principle, more reliable) results (the second column in two tables) of the recon- 
structed WIMP mass give nix ~ 120 GeV with a rough la upper (lower) bound 
of ~ 190 (80) GeV, or, equivalently, 

120t^° GeV . (2.4) 

Moreover, the combined results with two different form factors show not only a 
large overlap between ~ 85 GeV and 180 GeV, but also a good coincidence: 
comparing to the ^^40 GeV la statistical uncertainty and the ^^35 GeV overlap, 
the difference between two median values is ;< 10 GeV! This indicates that, for the 
first approximation of giving/ constraining the most plausible range of the WIMP 
mass, the uncertainty on the nuclear form factor could be safely neglected. 

2.2 Spin-independent WIMP-nucleon coupling |/pp 

Following the WIMP mass determination, I consider now the reconstruction of the 
SI WIMP couphng on nucleons l/pp (Shan 2011) with a '^*^Ge target^. 

In Figs. 2 I show the reconstructed squared SI WIMP-nucleon couplings and 
the lower and upper bounds of their la statistical uncertainties estimated by 
Eqs. (17) and (18) of Shan (2011) with an assumed (lOOilO GeV, labeled with 
the subscript "input") and the reconstructed (from Sec. 2.1, labeled with "re- 
con" ) WIMP masses. The commonly used value of the local Dark Matter density 
po = 0.3 GeV/cm^ and a larger value of po = 0.4 GeV/cm ' (Catena & UUio 2010; 
Salucci et al. 2010; Pato et al. 2010) as well as the elastic nuclear form factors 
given in Eqs. (2.1) and (2.2) have been used for estimating |/pp in the upper and 
lower frames, respectively. 

Among these results, the mean value and the overlap of two most plausible re- 
sults (estimated by using the reconstructed WIMP mass) give roughly (and some- 
how naively) a la range of 

l/pP - 9.0012-10 X 10-i« GeV-4, (2.5) 

or, equivalently, 

l/pl - 3.00;0:35 X 10-« GeV-2. (2.6) 

Since the reconstructed WIMP mass given iu Sec. 2.1 is rn,^ ^120 GeV, one can 
simply use the proton mass m-p to approximate the WIMP -proton reduced mass 



^Remind that the theoretical prediction by most supersymmetric models that the SI scaler 
WIMP couplings on protons and on neutrons are (approximately) equal: /p ~ /n has been 
adopted in the AMIDAS pa<;kage. 
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Results 

Reconstructed SI WEMP-nucleon coupling 
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Results 

Reconstructed SI WTMP-nucIeon coupling 
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Tlie reconstnicied squared spin-ind-ependent WIMP-nucleon coupling \fp\ arnS the lower and upper bounds of its 
lo statistical uncertainty estimated by Eqs. (17) and (18) of Ref. [3] with the input WIMP mass (with m overall 
uncertainty of 1%) and the reconstructed WIMP mass (estimated by the algorithmic; process described in Ref. [2]). 



Fig. 2. Tlie squared SI WIMP-nucleon coupling ]/p|'^ reconstructed with a ''^^Ge target. 
The commonly used value of the local Dark Matter density po = 0.3 GeV/cm"^ and a 
larger value of po = 0.4 GeV/cm"^ as well as the elastic nuclear form factors given in 
Eqs. (2.1) and (2.2) have been used for estimating |/pp in the upper and lower frames, 
respectively. 



771]., p and give a reconstructed SI WIMP-nucleon cross section as^ 

TT 



7SI 



^r,pl/pl 



"ip l/p 



4.31 



+1.01 

-0.69 



X 10"^ pb 



(2.11) 



Note that, since the expression for estimating l/pp (Eq. (17) of Shan (2011)) is a function 
of the (reconstructed) WIMP mass, for light WIMP mass, one has to use 



7SI 

XP 



1 

PO 



1 



1 



ySA^^ffn^ 



2Qif„r-(Q^i,) 



■lo 



(nix + "In) 



(2.7) 



where A is the atomic mass number of the target nucleus, £ is the experimental exposure. Then 
one has (cf Eq. (18) of Shan (2011)) 

- (4p) = 41 I , [1 + A(™x)] ' + A/-!-' (l/M.) 



+ '-^-"°"('"^'^/^-)ri + A(m,)1 
(rrix + "In) l J 



1/2 



Here I have used (Drees & Shan 2007) 

A/m = 

and defined 



FHQmin) 



A(m,) = 2 



(2.8) 

(2.9) 
(2.10) 
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Fig. 3. The reconstructed ratio between two SD WIMP-nucleon couplings, an/flp. As 
usual, the elastic nuclear form factors given in Eqs. (2.1) and (2.2) have been used for 
determining a-n/a^ in the upper and lower frames, respectively. 



2.3 Ratio of two spin-dependent WIMP-nucleon couplings a„/ap 

In Figs. 3 I show the reconstructed Qn/ap ratios and the lower and upper bounds of 
their Icr statistical uncertainties estimated by Eqs. (2.7) and (2.12) of Shan (2011) 
with n = 1 as well as by Eqs. (3.16) and (3.20) of Shan (2011) at the shifted 
energy points (Drees & Shan 2007; Shan 2011). A combination of ^^F + i27j 
targets has been used for the reconstruction of a^/up under the assumption that 
the SD WIMP-nucleus interaction dominates over the SI one (labeled with the 
superscript "SD"), whereas a third target of ^*Si has been combined with ^^F and 
-"^^^I for the case of the general combination of both SI and SD WIMP interactions 
(labeled with the superscript "SI + SD"). 

It can be found that, firstly, the "+ (plus)" solutions of the an/flp ratios given 
here are obviously too large to be the reasonable choice for a^/ap and the "— (mi- 
nus)" solutions should be the correct ones^. Secondly, although the reconstructed 
result under the assumption of the SD dominant WIMP interaction is in general 
larger than the (in principle more plausible) result obtained without such a prior 



Definitions and estimations of r(Qniin) and /„ can be found in e.g., Drees & Shan (2007, 2008). 

^Remind that, as discussed in Shan (2011), the correct choice from the "+" and "— " solutions 
can be decided directly by the values of the group spins of protons and neutrons of the used target 
nuclei, (^(p^n)}. 



Analyzing direct DM detection data by the AMIDAS website 



7 



Results ^^^^^^^^^^^^^^^^^^^^^^^^^^^^M 

Reconstructed ratio between the SD and SI WIMF-protoo cross sections 
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Tlie reoonslnjcted / "jtP ™tios and the lower and upper bounds of theii la statistical uncertaicties estimated 

byEqs. (50) aod (61) (wit n„/ Op estimated byEq. (57)) as weUas by Eqs. (55) and (69)of Ref. [4] at the shifted 
energy points (given by Eq. (38) of Ref. [4]). 

Reconstructed ratio between the SD and SI WEVIF-neutron cross sections 
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Results 








Reconstructed ratio between the SD and SI WIMP-proton cross sections ^^^^^^H 
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The reconslructed o^^^ 1 ratios and fhe lower and upper boiinds of their Iff statistical uncertaintie 
H by Eqs. (50) and (61) (with «□/ Op estimated by Eq. (57)) as well as by Eqs. (65) and (69) of Ref [4] at 
energy points (given by Eq. (38) of Ref [4]). 
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Fig. 4. The reconstructed ratios between the SD and SI WIMP-nucleon couplings, 
<T^p jj)/(j^p. As usual, the elastic nuclear form factors given in Eqs. (2.1) and (2.2) 



have been used for determining cr^J^ n)/''"xp upper and lower frames, respectively. 



assumption^ , one could still use the mean value and the overlap of these two results 
to roughly (and somehow naively) give a la range of 



--0.891°:^^ (2.12) 



2.4 Ratios of the SD and SI WIMP-nucleon couplings ct^d ^^/cr~^p 

In Figs. 4 I show the reconstructed '^^p n)/'^xp ^'^tios and the lower and upper 
bounds of their la statistical uncertainties estimated by Eqs. (3.9), (3.10) and 



^See also discussions in Sec. 2.4. 
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(3.21) of Shan (2011) (with a^/a^ estimated by Eq. (3.16) of Shan (2011)) as well 
as by Eqs. (3.25) and (3.29) of Shan (2011) at the shifted energy points. 

By using the data sets of ^^F, ^^''I and ^^Si targets (labeled with the superscript 
"XYZ") or combining that of ^^Na or ^^^Xe with the (common) data set of ^^Ge 
(labeled with the superscript "XY"), one can use the mean value and the overlap 
of these two results to roughly (and somehow naively) give a Itr range of 



SD 

-^^9.6ll^:ixl0^ 



.SD 
X" 
-SI 
XP 



5.45t^;^^ X 10^ . 



Then, firstly, from these results one can further obtain that^' ^ 



(2.13) 



0.75 



+0.15 
-0.23 ■ 



(2.16) 



Secondly, combining the results in Eq. (2.13) with tr^p given in Eq. (2.11), one can 



also obtain that^ 

SD 



XP - 4.14tl:^^ X 10-3 



(JC;C =i X iU " pb , 

These results give in turn that^ 



U.iU8_o.o20 I 
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^Here I have used 
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(2.14) 



(2.15) 



since the 1<t uncertainties given in Eq. (2.13) are not the exact but only rough estimates from 
the overlaps of two results given in Figs. 4. 

'^Remind that the results given in the second and third columns of the tables in Figs. 4 axe 
reconstructed with the On/op ratio given in the last columns of the tables in Figs. 3. 

^Here I have used 



'"K(p,n)j - [Kp] K(p,n)/<^xpj + K(p,n)/<^: 

and neglected the correlation term in the bracket: 
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by assuming that two independent data sets with the Ge target and other two independent 

XP 



data sets with the ^*Si target have been used for determining (T^p and n)/"^!?" 



Since 



1 '2^F"^?,(p,n)l»(p,n)l^ 



(2.20) 
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On the other hand, one can also use the reconstructed On/op ratio given in 
Eq. (2.12) and one of the two results given in Eq. (2.19) to obtain that^° 

^ 2.9712-05 X 10-3 pb , ~ 3.28t2;24 ^ IQ-^ pb . (2.26) 

These results can also give that 

|ap| ~ o.o9itn!o , l«n| ^ o.ogetHs? • (2-27) 

It can be found that, not surprisingly, the statistical uncertainties on the re- 
constructed cr^l^ given in Eq. (2.26) are 2 or 3 times larger than those given 
in Eq. (2.19): Since '''^(p n) /"''xp reconstructed with the F + I + Si combination 
involve already the reconstructed a^/a^ ratio given in Eq. (2.12), the uncertain- 
ties on cr^^ j^-j given in Eq. (2.26) are thus overestimated. Secondly, although the 
reconstructed cr^p and given in Eqs. (2.19) and (2.26) have overlaps, these 
results seem not to match to each other very well; given in Eq. (2.26) is even 
larger than there although the an/flp ratios given in Eqs. (2.12) and (2.16) are 
< 1. One possible explanation is that the On/cip ratio given in Eq. (2.12) would be 
overestimated. This can be seen by comparing the fln/ftp ratio given in Eq. (2.12) 
to that given in Eq. (2.16) estimated (somehow independently) by the results given 
in Eq. (2.13). 

Nevertheless, the analyses given here show that, firstly, once one can estimate 
the SI WIMP-nucleon coupling/cross section, |/p| or a^p, and (one of) the ratios 
between the SD and SI WIMP-nucleon cross sections, and/or the ratio between 
two SD WIMP-nucleon couplings, the other couplings/cross sections could in prin- 
ciple be estimated. Secondly, although the method under the assumption of the 



and ~ 120 GeV, one has 
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^"Here I have used 
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^K(p,n)j -<^x(p.n) 



1/2 



(2.21) 



(2.22) 



(-^Fn,P)) /«(n.p))' +4.^(«n/ap) /(a„/ap)^ 
and neglected the correlation term in the bra<;ket: 

=F4cov(aSD _pj,a„/ap) /(<^xfn,p)) (»n/ap) , 

since the la uncertainties given in Eq. (2.12) as well as in Eq. (2.19) are not the exact but only 
rough estimates from the overlaps of the results given in Figs. 2, 3 and 4. The "— (-1-)" sign in 
Eq. (2.25) is for the case with protons (neutrons). 



(2.24) 



(2.25) 
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SD dominant WIMP interaction would overestimate (or underestimate, depending 
on the combination of the used targets (Shan 2011)) the a^/up ratio, the recon- 
structed result (s) could still be useful for at least determining the correct sign of 
Qn/ttp. Moreover, the WIMP couplings/cross sections estimated in different way 
would be self-cross-checks to each other and the (in) compatibility between the 
reconstructed results would also help us to check the usefulness of the analyzed 
data sets offered from different experiments with different detector materials. 

3 Input setup for generating pseudodata 

In Table 1 I give finally the input setup for generating the pseudodata sets used 
in the analyses demonstrated in the previous section. For comparison, the recon- 
structed results shown in the previous section are also summarized here. 

It can be found that, firstly, not only the WIMP mass given in Eq. (2.4) and 
the result reconstructed with the input nuclear form factor (lower frame of Figs. 1), 
but even the mass reconstructed with the "wrong" form factor (upper frame) can 
match the input WIMP mass very well: the deviations between the input and the 
reconstructed values are only ~ 13% (with the wrong nuclear form factor) or even 
only ^ 6% (with the input one). As discussed earlier, this indicates that, for the 
first approximation of giving/constraining the most plausible range of the WIMP 
mass, the uncertainty on the nuclear form factor could be safely neglected. 

Secondly, all WIMP nucleon couplings/cross sections as well as the ratios be- 
tween them have also been reconstructed with only ~ 5% to ^ 40% deviations 
from the input/theoretically estimated values. Although the SI WIMP coupling 
I /pi estimated with the input (larger) local Dark Matter density (lower frame of 
Figs. 2) is underestimated (Shan 2011), one can at least give an upper bound on 
|/p|. Meanwhile, although the Un/ap ratio given in Eq. (2.12) is overestimated, in 
Sec. 2.4 w(; have demonstrated that by combining different methods for estimating 
different (ratios between the) WIMP couplings/cross sections, one could in prin- 
ciple observe/confirm the (in)compatibility between these results and probably 
correct the reconstructed values. 

Moreover, for generating pseudodata, we have used the shifted Maxwellian 
velocity distribution: 



with the Sun's Galactic orbital velocity vq = 230 km/s; v^. is the time-dependent 
Earth's velocity in the Galactic frame: 



the date on which the Earth's velocity relative to the WIMP halo is maximal has 
been set as tp = 140 d. Although these values for the astronomical setup are 
non-standard, we would like to stress that, firstly, for using the AMIDAS package 




(3.1) 




(3.2) 
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Property 


Reconstructed value 


Input/Estimated value 


Remarks 




UOtll GeV 


130 GeV 




XP 


4.3ltJ:°i X 10-9 pb 


4 X 10-9 pb 


/n /p 


l/pp 


9mtl:ll X 10-18 GeV-^ 


9.305 X 10-18 GeV-^ 


t 


I/pI 


Smto ll X 10-9 GeV-^ 


3.050 X 10-9 GeV-^ 


t 


Op 




0.1 




On 


omits 


0.07 






0.89t°i^, 0.75t°:^^ 


0.7 




SD 
XP 


4.14+}-4^ X 10-3 pb 


3.51 X 10-3 pb 


t 








+ 

T 


SD / SI 
XP ' XP 


9-61±3:28 X 


8.77 X 10^ 


t 


^SD/ SI 
Xn / XP 


5-45l2:75 X 10^ 


4.30 X 10^ 


t 


^^li(Q) 




i^li(Q) in Eq. (2.2) 








F|s(Q) inEq. (2.3) 




Po 




0.4 GeV/cm^ 




tp 




140 d 




^cxpt 




300 d 




Wo 




230 km/s 




^max 




600 km/s 




(^expt ) 




226.6 km/s 





Table 1. The input setup for generating the pseudodata sets used in the analyses demon- 
strated in this article. The theoretically estimated values and the reconstructed results 
shown in the previous section are also given, f: estimated for 130 GeV m^- 



and website to analyze (real) data sets, one needs only the form factors for SI 
and/or SD WIMP-nucleaus cross sections, prior knowledge/assumptions about 
the WIMP velocity distribution fi{v) and local density po (except the estimation 
of the SI WIMP-nucleon coupling |/pp) are not required. Secondly, as shown in 
the previous section, such non-standard values would not affect the reconstructed 
results. 

4 Summary 

In this article I demonstrated the data analysis procedures for extrating WIMP 
properties by using theoretically generated pseudodata for different target nu- 
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clei. As an extension as well as the complementarity of our earlier theoretical 
works, I combined reconstructed results of the (ratios between different) WIMP 
couplings/cross sections on nucleons to estimate each individual coupling/cross 
section. Hopefully, the AMIDAS package and website as well as this demonstration 
can help our experimental colleagues to analyze their real direct detection data 
in the near future and to determine (at least rough ranges of) properties of halo 
Dark Matter particles. 
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